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Abstract 

A Monte Carlo sensitivity analysis is used to build a plasma chemistry model for octacyclofluorobutane (c- 
QF 8 ) which is commonly used in dielectric etch. Experimental data are used both quantitatively and 
qualitatively to analyze the gas phase and gas surface reactions for neutral radical chemistry. The 
sensitivity data of the resulting model identifies a few critical gas phase and surface aided reactions that 
account for most of the uncertainty in the CF 2 and CF radical densities. Electron impact dissociation of 
small radicals (CF 2 and CF) and their surface recombination reactions are found to be the rate-limiting 
steps in the neutral radical chemistry. The relative rates for these electron impact dissociation and surface 
recombination reactions are also suggested. The resulting mechanism is able to explain the measurements 
of CF 2 and CF densities available in the literature and also their hollow spatial density profiles. 

I. Introduction 

Octacyclofluorobutane (c-C 4 F 8 ) is a commonly used dielectric etch gas with rich polymerizing properties. 
The differences in deposited C x F y film characteristics and their deposition rates on different materials 
enables the etch process to proceed at disparate rates offering the much needed etch selectivity. In order to 
meet the stringent requirements of the next generation semiconductor integration techniques, it is critical 
that we gain enough quantitative insight into these processes to be able to build models with sufficient 
fidelity. It is believed that the precursors of the polymer film formation are the neutral radicals produced in 
the plasma via electron impact dissociation of feedstock gas molecules and their fragments. It has also 
been suggested that low energy ion bombardment aids in the formation of polymer film whereas high 
energy ion bombardment at a biased electrode adversely affects the film formation due to sputtering. 1 
Therefore, a better quantitative understanding of the plasma chemistry that forms these radicals and ions is 
an essential step toward understanding these surface mechanisms. 

Due to the large size of the c-C 4 F 8 molecule, it forms a wide variety of ions and neutrals due to electron 
impact ionization and dissociation. 2 These fragments then undergo further ionization and dissociation into 
smaller ions and neutral radicals due to electron impact. Negative ions are also formed due to attachment 
and dissociative attachment reactions. As ions and neutrals collide, they alter the relative populations of 
ions by charge exchange. In addition, all species interact with the chamber walls. Most ions are neutralized 
at the surfaces. However, a neutral radical may also be adsorbed on a surface and recombine or react with 
other surface species and may be released as a radical or a stable species. The rates at which these gas- 
phase and gas-surface reactions occur are very difficult to measure. Currently, only the rates of electron 
inpact ionization reactions of stable species ( e.g . : c-QFg, CF 4 , CHF 3 , etc.) are known with good certainty. 2 * 
4 Electron inpact ionization rates of only a handful of neutral radicals (CF 3 , CF 2 , and CF) have ever been 
measured. 5 Among the gas-phase reactions, the biggest roadblock toward building a quantitative model is 
the lack of credible measurements of electron impact dissociation of stable species and radicals into smaller 
neutral radicals. This is difficult to measure due to difficulties in tracking neutral radical products. As a 


result, neutral dissociations reaction rates are scarce in the literature and it is not uncommon for the rates 
obtained by different groups to differ from each other by orders of magnitude. 

Clearly, our chemical kinetics database needs are vast and it is unlikely to be fulfilled in the near future. 
Therefore, we must rely on alternative means to reduce the uncertainties in our model calculations. One 
such way is to perform a sensitive analysis of the plasma system and identify a few key processes that are 
responsible for the majority of the uncertainty in the model. Once these major sources of uncertainty are 
identified, comparisons with a variety of plasma diagnostic data can help reduce some of the uncertainty. 
Also, an uncertainty analysis identifies critical rate-limiting and dominant reaction mechanisms which can 
be the focus of future investigations, essentially narrowing the database needs. In the past, we performed 
uncertainty analysis for a CHF 3 plasma system. 6 

In this paper we analyze c-C 4 F 8 chemistry in a high density plasma reactor. There have been a few 
multidimensional reactor modeling studies of c-C 4 F 8 . 7-9 Due to its extensive use in dielectric etch, some 
experimental diagnostics data are also available in the literature. Hebner and Abraham 10 and Hebner 11 
present plasma densities and some neutral radical densities in an inductively coupled Gaseous Electronics 
Conference (GEC) cell. They also presented some radial profiles of CF 2 radicals, the dominant neutral 
species in the high density chamber. Waters 12 also measures CF 2 and CF radicals and electron densities in a 
GEC cell during oxide and photoresist etch. There have also been some measurements of neutral radical 
profiles in ECR discharges. 13,14 Almost all neutral density profiles measured by variety of experimental 
groups show a hollow profile for CF 2 density, i.e. CF 2 is larger near the chamber walls as compared to the 
center of the discharge. This essentially means that there is a net production of CF 2 at the chamber walls 
and a net destruction in the bulk plasma. This has also been observed in plasma beam measurements for 
other fluorocarbon systems. 15,16 CF X profile measurements made in low density capacitively coupled CF 4 
discharges also show a net production of CF 2 radicals at the electrodes. 17,18 In this work, we will show that 
by comprehensively comparing our data to these experimental measurements and with the help of 
sensitivity analysis we can build a CF 2 formation and consumption model by isolating the rate-limiting and 
dominant steps. The relative magnitudes of the rates of these steps, if not the rates themselves, may be 
deduced from this analysis. 

The rest of the paper is organized as follows. In Sec. II we discuss the c-C 4 F 8 plasma chemistry and the 
sources of the rate constants. Section III outlines the details of the plasma simulation and comparisons with 
experimental data. The sensitivity analysis technique is discussed in Sec. IV and in Sec. V we discuss the 
results of the sensitivity analysis and discuss their implications. Finally we make some concluding remarks 
in Sec. VI. 

n. c-C 4 F 8 Plasma Chemistry Model 

Christophorou and Olthoff 2 present a comprehensive review of the electron collision data for c-C 4 F 8 . In this 
section we present a discussion of the available data that is used in this work. Also, since in high density 
plasmas most of the feedstock c-C 4 F 8 is dissociated into smaller radicals, the electron collisions with these 
smaller radicals are more important. 

Electron impact dissociative ionization: 

The electron impact of c-C 4 F 8 produces a variety of ions; C 3 F 5 + , C 2 F 4 + , CF 3 + , CF 2 + , CF + , etc. There have 
been a few measurements of electron impact ionization cross sections. 19,20 In this work we use the data 
presented by Jiao et al. 19 The reactions considered are listed in Table 1. Although the data obtained by Jiao 
et al. are not in very good agreement with the data obtained by Toyoda et al., 20 these reactions are not 
expected to strongly influence the neutral radical composition of the plasma as long as the total ionization 
cross sections are consistent. Also, most of the ionization in a high density plasma occurs via electron 
impact of dominant neutral radicals, in this case CF 2 and CF, for which we use the fast neutral beam data of 
Tamovsky et al. 5 We also include atomic ionization reactions as shown in Table 1. 

Electron impact dissociation into neutrals: 


\ 



It is generally believed that since c-QFg is a dimer of C 2 F 4 , it primarily dissociates into two C 2 F 4 
molecules. Bauschlicher and Ricca 21 using ab initio bond energy computations have confirmed that c-C 4 F 8 
is much more likely to dissociate by breaking two opposite C-C bonds forming C 2 F 4 (1.95 eV) rather than 
two adjacent C-C bonds (3.26 eV) and forming C 3 F 6 +CF 2 . The breaking of a C-F bond and forming QF 7 +F 
is even more unfavorable (4.7 eV). This is has been experimentally seen in shock tubes where thermal 
dissociation of c-C 4 F 8 mostly produces C 2 F 4 with only small amounts of C 3 F 6 . 22,23 Also, the activation 
energy of dissociation computed from the ab initio study agrees well with the experimental estimate, 23 
lending further credibility to the computations. Bauschlicher and Ricca 24 further show via ab initio heat of 
formation computations that for C x F y (x>l) radicals, pulling a CF 2 group out is easier than other 
dissociation pathways. This implies that C 2 F 4 and C$F 6 will also most likely dissociate by releasing a CF 2 
group. Therefore CF 2 is readily formed in a c-C 4 F 8 plasma. The CF 2 radicals are also further dissociated 
into CF, C, and F radicals by electron impact. The reactions considered in this work are shown in Table 1. 
A few authors have estimated the total dissociation rates of c-C 4 F 8 by electron impact at specific electron 
temperatures. 25,26 In this work we use the Arrhenius expression of Font et al. 9 which is estimated from 
swarm measurements. We divide the total dissociation rate into two major dissociation pathways with a 9: 1 
branching ratio. 

e+C 4 F 8 ->2C 2 F 4 +e 
e+C 4 Fg — > C3F0 +CF 2 +e 

This ratio is chosen based on mole fraction ratios of C 2 F 4 and C^ seen in thermal decomposition of c- 
QFg. 23 As we will find later, the neutral radical densities (CF 2 , CF, etc.) are not very sensitive to this 
choice. The dissociation rate of C 2 F 4 is taken from Font et.al. 9 as well. However, this rate has not been 
estimated experimentally and is likely to be somewhat uncertain. The dissociation rates of remaining 
radicals haven been estimated from other systems, as tabulated in Table 1. More on dissociation cross 
sections will be discussed later. 

Ion-molecule charge transfer: 

Ion-molecule charge transfer reactions alter the relative densities of various ions. However, since the 
densities of dominant radicals are orders of magnitude larger than that of the ions, they are not significantly 
affected by these reactions. We include a variety of charge transfer reactions outlined by Font et.al. 9 The 
rates of these reactions are mostly determined by the Langenvin formula. The uncertainties in these rates 
will directly affect the relative densities of various ions, whereas the neutral radicals which are the focus of 
this article are unlikely to be affected very much. In very high density plasmas ( n e > 10 12 cm' 3 ), the ion 
composition of the plasma can affect the neutral composition as the wall recombination of ions becomes a 
significant contributor to the neutral species mass balance. This is unlikely to be the case in the parameter 
space that we explore and will also be evident from the sensitivity data. 

Attachment , detachment , and ion-ion recombination: 

At low electron energies ( < 1 eV) the c-C 4 F 8 - negative ion has been observed in electron beam studies. 27 In 
a strongly dissociated high density discharge negative ions are most likely to form via dissociative 
attachment of CF X radicals. The rates of these dissociative attachment reactions tend to be uncertain since 
there are no direct measurements of these individual reactions. Attachment reactions have indirect 
influence on neutral radical densities via changes in plasma properties, such as plasma potential, density 
and temperature. In this work we adjust the rates of CF X dissociative attachment reactions to get good 
agreement with experimental data presented in the next section. A few ion-ion recombination reactions 
listed by Font et al. 9 have also been included. Their influence on the neutral radical densities is expected to 
be minimal. 

Gas-surface reactions: 

The current state of knowledge of gas-surface reactions in fluorocarbon plasmas is much more limited than 
the gas-phase reactions. In most fluorocarbon plasmas, etching and deposition processes proceed 
simultaneously via a variety of surface reaction mechanisms, such as physical and chemical absorption, 
surface diffusion, desorption, ion impact processes. A complete modeling based on first principles is 



currently not feasible due to the lack of understanding of the fundamental processes. Therefore, in this 
paper we rely on phenomenological modeling of gas surface reactions. 

In many fluorocarbon plasmas radical density profiles have been measured using LIF and UV 
absorption. 11 ’ 13,14,17 Based on these measured profile gradients we can assess whether chamber walls act as a 
source or a sink for a particular neutral species. For CF 2 radicals, a dominant species in c -C 4 F 8 plasmas, the 
density profiles shows a hollow shape, i.e. a larger density at the chamber walls as compared to the one in 
the bulk plasma. A part of this density gradient may also be due to the gas temperature profile, but it is 
unlikely that a gas temperature gradient can be large enough to cause such large CF 2 density gradients as 
were measured. Therefore, it was reasonably concluded by many authors that there is a net production of 
CF 2 at the chamber walls. For CF, however, there appears to be a density gradient directed away from the 
chamber walls, indicating a possible consumption at the walls. The measured CF density profiles of 
Nakamura et al. 13 and Suzuki et al. 14 also show a hollow profile, but unlike the CF 2 density, the CF density 
gradient changes direction near the chamber walls. It remains unclear what mechanisms are responsible for 
CF 2 production at the chamber walls. It has been speculated that larger polymeric neutrals may deposit on 
the chamber walls forming the amorphous film which are then desorbed as CF 2 . 18 However, later 
measurements have shown that the density of large polymeric neutrals is too low in high density plasmas to 
justify a large CF 2 flux from the walls. 28 It has also been suggested by observations in plasma beam studies 
that CF 2 radicals may be produced at the walls by CF+F recombination. 15,16 Based on the above 
observations, we believe that it is highly likely that CF radicals recombine on the fluorinated surfaces of the 
chamber forming CF 2 which is readily desorbed. The film formation occurs at a much slower time scale 
than the production of CF 2 . Therefore, it is not unreasonable to assume that the slow process of sticking of 
radicals leading to the fluorocarbon film formation is only a subset of the variety of dissociation, exchange, 
and recombination processes that may be occurring on these surfaces. The details of these processes is 
obviously unknown, but in this work we start by assuming that following processes occur at all surfaces 

cf + f->cf 2 (2) 

O + F — ^ OF 

It is possible that further recombination reactions forming CF 3 and CF 4 also occur at these surfaces, but 
since these species are not found with appreciable densities in c-C 4 F 8 plasmas, we will ignore this 
possibility. The rates of surface aided recombinations shown in (2) are obviously unknown, but in this work 
we attempt to estimate them via sensitivity analysis and comparisons with experiments. 

In addition to the above recombination reactions, we assume that all ions recombine at the walls with a unit 
probability. The complete list of all gas-phase and surface aided reactions are shown in Table 1. 


HI. Volumetric Plasma Simulations and Comparisons with Experiments 

In this section we briefly describe the plasma reactor scale simulations. We use the SAMPR code 32 to solve 
volume averaged conservation equations for each neutral and ionic species with surface reaction boundary 
conditions as shown in (2) along with ion neutralization reactions. In addition we solve for electron and gas 
temperatures. All ions are assumed to exit the bulk with individual Bohm speeds. The reactions shown in 
Table 1 are included. In order to establish a baseline plasma chemistry model we begin by comparing and 
calibrating our model results with two sets of experimental measurements of plasma and neutral radical 
densities in pure c-C 4 F 8 plasma in an inductively coupled GEC cell. The first set of data we consider is by 
Waters 12 who used infrared diode. laser absorption spectroscopy to measure CF 2 and CF densities at a 
variety of pressures, source powers and bias powers. The data were condensed by Waters in the form of ^ 
fitted expressions which are used in this paper to validate the results of our parametric study. Waters also 
presented Langmuir probe measurements of electron density and plasma potential. The second set of data is 
by Hebner 11 and Hebner and Abraham 10 who presented CF 2 density on a relative scale using LIF in a GEC 
cell. They also presented electron and F- ion densities. 


Comparisons with the experiments of Waters: 



We compare the CF 2 and CF density variations with electron density at pressures of 8,10, and 20 mTorr in 
Fig. 1. The flow rate is fixed in this experiment at 26.1 seem of pure c-C 4 F 8 . We plot radical densities 
versus electron density, instead of source power. This is done to eliminate the uncertainties that cause 
errors in source power and electron density relationship. These uncertainties may be due to the lack of 
knowledge of total plasma power, other power consumption mechanisms such as vibrational and electronic 
excitations, and power losses in the RF circuitry. Due to difficulties involved in modeling these 
mechanisms, electron densities are often obtained in models by assuming that only a certain fraction of the 
source power is actually used to sustain the plasma. In electronegative plasmas the prediction of electron 
density is further complicated by attachment reactions whose rates are usually very uncertain. In order to 
eliminate these uncertainties and focus only on processes that generate (or consume) neutral radicals we 
plot radical densities versus electron density, thereby ensuring that plasma conditions are identical in 
experiments and simulations. 





Figure 1. Comparison of the model CF 2 and CF densities with the fits of the experimental measurements by 
Waters. 12 (a) 8 mTorr, (b) 10 mTorr, and (c) 20 mTorr. 

In Fig. 1 the experimental data is plotted using the fitted curve obtained by Waters 12 and the model results 
obtained in this work by symbols. For all three pressures considered CF 2 density falls as the plasma density 
rises (with increasing source power). This is primarily due to increased dissociation of CF 2 in a high 
density discharge. In low density capacitive discharges often CF 2 density increases due to an increased 
dissociation of the feedstock gas and an increased production of CF 2 at the powered electrode. In the 
inductively coupled discharge with tens of milliseconds of residence time, almost all of the feedstock gas is 
already dissociated into smaller radicals, such as CF 2 and CF. A further increase in plasma density causes 
the CF 2 molecules to dissociate into CF and atomic species resulting in a rise of CF density (see Fig. 1). 
The models results agree very well with the experiments at 8 and 10 mTorr. The slopes of the radical 
densities match particularly well with the experiments. However, as the pressure rises from 8 to 20 mTorr 
the model predictions of CF 2 and CF densities become larger than the ones measured in the experiments. 
This is very likely because of the spatial variations in the density profiles. As discussed before CF 2 and CF 
density profiles exhibit a hollow profile with low density in the middle of the discharge compared to the 
periphery. Therefore, it is expected that the experimental densities which are line averaged will be smaller 
than the model predictions that are volume averaged. As the pressure increases, the density variations (the 
hollowness) become more pronounced resulting in a larger difference between model and experiments. 






Figure 2 . Comparison of the model (a) CF 2 density and (b) F- ion density with the values measured by 
Hebner ll and Hebner and Abraham . 10 

Comparisons with the experiments of Hebner and Abraham: 

Hebner 11 gives CF 2 density measurements on a relative scale in C-C 4 F 8 plasma. The data for varying power 
is given only for 10 mTorr. In order to plot the CF 2 density against the electron density we use 
corresponding values given by Hebner and Abraham . 10 Our model data agrees very well with the CF 2 
density variation with increasing electron density within the experimental scatter as shown in Fig. 2(a). We 
normalize the experimental CF 2 data with the model CF 2 density at the lowest electron density point. 
Again, the slope of the CF 2 density, which is the only relevant quantity in this case, is predicted by the 
model well. 

Hebner and Abraham 10 also presented F‘ negative ion density data which is found to be the dominant 
negative ion. We find that F" is critically dependent on the dissociative attachment reactions involving CF 2 
and CF radicals. The rates of these reactions are adjusted to attain the best possible match with the 
experiments as shown in Fig. 2(b). 

IV. Sensitivity Analysis Technique 

The goal of this sensitivity analysis is to isolate the most sensitive gas phase and surface reactions that 
determine the CF 2 and CF radical densities in a high density c-C 4 F 8 discharge. Since we obtained very good 
agreement with the available experimental data in the previous section using our baseline model, the 
sensitivity analysis will allow us to place confidence levels on some of the critical reaction rates. Ideally, 
this methodology could be used to reduce uncertainties in our reactions rates or at least establish relative 
rates of certain critical reactions. 

Since a plasma chemical system is inherently non-linear, we use a Monte Carlo technique where the 
reaction rate for each reaction is varied simultaneously and independently based on a probability 
distribution function centered at the baseline rate coefficient or the reaction probability for that reaction. 
The baseline rate coefficients and probabilities are shown in Table 1. A Gaussian probability distribution 
of Iogio£ is chosen for each reaction with a half width of half-order of magnitude (Alogio k = 0.5). The 
width of the Gaussian is chosen based on the uncertainties of most of the reaction rates. An infinitesimally 
small Gaussian width would reduce the technique to a local sensitivity analysis, which is frequently 
undertaken. However, the large variability in reaction rate coefficients makes the non-linear effects more 
important. 



The steps involved in this methodology are as follows. (1) Random numbers are used to generate rate 
coefficients based on the probability distribution. (2) The plasma model SAMPR is run using these rate 
coefficients and species, and plasma densities are collected. (3) In this work steps (1) and (2) are repeated 
tens of thousands of times to gather a large enough sample to drawn meaningful correlations. (4) 
Correlation coefficients are obtained for each input-output pair. (5) Rate limiting steps are identified. For 
details on these steps the reader is referred to Ref. 33 and our earlier work on CHF 3 sensitivity analysis. 6 

V. Sensitivity Results and Discussion 

In this section we discuss the sensitivity analysis data obtained for pure c-C 4 F g plasma in a GEC cell at 10 
mTorr with 26.1 seem of flow rate and 500 W of plasma power. The bias power is set to zero. The data are 
obtained in the form of scatter correlation plots. In this work we ran 10,000 runs which are found to be 
sufficient to describe the statistics. We use linear regression analysis to obtain correlation coefficients and 
sensitivity rankings. 

CF 2 density: 

For CF 2 density we find that only a handful of reactions are important. Figure 3 shows the correlation 
coefficients of CF 2 density with all reaction rate coefficients. The reactions with large magnitudes of 
correlations coefficients are the most sensitive reactions. The results show that the two main producers of 
CF 2 are the surface recombination of CF+F and the electron impact dissociation of C 2 F 4 (which is produced 
from the feedstock gas) in the bulk. The primary consumers of CF 2 are the electron impact dissociation and 
ionization reactions in the bulk, which also have large correlation coefficients. These findings are 
consistent with hollow profiles of CF 2 density that are observed experimentally, 13,14,11 which result from a 
net consumption of CF 2 in the bulk and a net production at the surfaces. In the bulk plasma, although 
feedstock gas dissociation is a large producer of CF 2 (reactions 23 and 27 in Table 1), there is a net 
consumption of CF 2 in the bulk due to fast electron impact dissociation and ionization of CF 2 (reactions 17 
and 31 in Table 1). As can be seen in Fig. 3(a), CF 2 is in not sensitive to the rates of the remaining gas- 
phase reactions in the range over which the rates are varied (+ half order of magnitude). This implies that 
our current knowledge of these reactions rates is sufficient and further reduction of uncertainties in the rates 
of these remaining reactions does not make the model results more reliable. The rate of the surface 
recombination reaction CF+F-*CF 2 is, however, strongly correlated with the CF 2 density and is a major 
producer of this radical. Fig. 3(b) shows the correlation plot of the CF 2 density with the rate of CF+F 
surface recombination reaction. An order of magnitude change in the recombination rate changes the 
density of CF 2 by more than a factor of 2. 




CF+F Recomb. CoefF. 


Figure 3 (a) Sensitivity of CF 2 density on all gas phase and gas surface reactions considered in Table 1. (b) 
CF 2 density correlation with the coefficient of CF+F recombination at the surface. 



CF density: 

For CF density determination only two reactions are important. The production-consumption cycle of CF 
appears very simple. CF is produced in the bulk primarily by dissociation of CF 2 as seen Fig. 4(a). A small 
negative correlation with reaction 17 is because as more CF 2 is ionized, less CF 2 remains for eventual 
production of CF. The primary consumer of CF is the wall recombination of CF+F. The recombination of 
C+F at the wall is also considered, but since the density of C atoms is very low, this reaction appears 
inconsequential as far as CF density is concerned. Since CF is produced in the bulk and consumed at the 
surface, the CF density profile is not expected to be hollow, unless there is a peaked gas temperature 
profile. Reliable gas temperature measurements are not available in the literature. Rotational temperatures 
measurements that are often presented show significant peaked profiles, but the gas may not be in local 
rotational equilibrium with the translational temperature. 34 In any case, experiments show both peaked and 
hollow profiles for CF density, but the gradient of CF density near the wall in the measurements so far 
indicates a surface loss of CF. 13,14 Figure 4(b) shows the correlation plot of CF density with CF+F 
recombination at the wall. The correlation is clearly strongly negative and nonlinear, which indicates that 
the sensitivity of CF density to the recombination rate changes over the range of variability. 



Figure 4 (a) Sensitivity of CF density on all gas phase and gas surface reactions considered in Table 1. (b) 
CF density correlation with the coefficient of CF+F recombination at the surface. 



Flow Rate (seem) 



Figure 5. Effect of flow rate on the sensitivity of CF 2 density on major gas phase and gas surface 
production reactions. 

Rate-limiting steps: 

For CF X radical production we find that the feedstock gas dissociation, CF 2 dissociation and CF+F wall 
recombination are the rate-limiting steps. As discussed previously, c-C 4 F 8 dissociates via electron impact 
primarily into two C 2 F 4 molecules, which then dissociate into four CF 2 molecules. From the best rates 
available in the literature it appears that the first step in this series (e+c-C 4 F 8 ^2C 2 F 4 +e) is faster than the 
second step (e+C 2 F 4 — >2CF 2 +e), and is not a rate-limiting step. The second step is one of the rate-limiting 
steps for CF 2 production in the bulk, but only if the discharge is not flow limited. At low flow rates the 
feedstock gas rapidly converts into CF 2 (see Fig. 5) and C 2 F 4 dissociation ceases to influence the CF 2 
density reflected by its low correlation coefficient. In this case, the rate of CF 2 production in the bulk is 
controlled by the rate of gas injection. It is worth emphasizing here that despite CF 2 production in the bulk, 
there is a net consumption of CF 2 in the bulk due to fast CF 2 dissociation and ionization processes. The rate 
of CF 2 dissociation is currently unavailable in the literature. In this work we estimate this rate from typical 
dissociation rates of other molecules. Figure 6 shows measured electron impact dissociation rates of some 
diatomic and polyatomic molecules. In the temperature range of interest it appears that most dissociation 
rates safely he within an order of magnitude from each other. We surmise that CF 2 dissociation rate also 
probably falls within this one order of magnitude band of uncertainty and assume that the CF 2 dissociation 
rate is same as that of CF 4 dissociation. Given this choice, the recombination coefficient of the other rate- 
limiting step, the CF+F— »CF 2 (surface), is set to 0.2 since it gave the best agreement with the experiments 
seen earlier in the paper. The validity of the rates chosen for CF 2 dissociation and CF+F recombination is 
currently unknown, but the relationship between these rates is probably correct since it is the balance 
between these two reactions that gives us CF 2 and CF densities in the discharge that are found to agree well 
with the experiments. The role of ion impact on the CF+F surface recombination rate is not explored in this 
work. Since all surfaces were unbiased, the effect of ion impact on surface recombination rates may be 
minimal or constant over the rage of conditions explored in this work. 



Figure 6. Electron impact rate coefficients of a few diatomic and polyatomic molecules. The rates for CF4, 
CHF3, CO, and 02 dissociation reactions are obtained from Refs. 35,6,36,37 respectively. 


VL Concluding Remarks 

From the Monte Carlo sensitivity analysis we find that CF 2 and CF densities in a high density plasma 
reactor are primarily determined by the balance between the rates of electron impact dissociation of CF X 
radicals and their surface recombination reactions. The rates of feedstock gas dissociation and neutral 
radical ionization play secondary roles in determining the radical densities. Since the electron impact 
dissociation of CF X and surface recombination are the rate limitin g steps, we estimate their relative rates by 




comparing the radical densities with the experimental data. Using the estimated rates and the rates obtained 
from the literature we are able to obtain excellent agreement with the radical density measurements of 
Waters 12 and Hebner. 11 The surface recombination reactions for neutral radicals are also able to explain the 
hollow profiles of CF 2 obtained by various authors. 11,13 ’ 14 The neutral chemistry model can be further 
refined if experimental measurements of spatial profiles of gas temperature and atomic species densities 
become available in the future. 
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Table 1. Gas phase reactions considered in the c-C 4 F 8 plasma model. The symbol y represents the reaction 
probability for the gas-surface reactions. 


Id 

Reaction 

A (cm 3 /s) 

T 1 

D (eV) 

AE 

(eV) 

Ref 

Electron impact excitation of c-C^g: 

1 e + C 4 F 8 -> C 4 F 8 (v) + e 3.390xl0' 8 

-1.093 

0.63 

0.12 

9 

Electron Impact Ionization: 

^ G + C 4 F 3 — ^ CF 3 + C 3 Fg^ + 2e 

6.655x10-® 

0.410 

18.71 

11.35 

9,19 

3 

e + C 4 F q — > ^ 2^4 + 2e 

5.568x10'® 

0.470 

17.46 

12.0 

9,19 

4 

e + C 4 F 8 — > C 2 F 5 + OF 3 + 2s 

2 . 688 x 10 '® 

0.379 

22.30 

12.58 

9,19 

5 

s + C 4 F 8 — > C^Fg + OF 2 + 2s 

4.840x10'® 

-0.026 

27.03 

14.63 

9,19 

6 

s + C 4 F 8 — > C^Fg + F + CF + + 2e 

1.279x10'® 

0.654 

20.61 

17.6 

19 

7 

e + C 3 F 6 — y F + G 3 F^ +2s 

3.025xl0' 9 

0.874 

16.41 

13.0 

From C2F4 

8 

e + CgFg — ► C 3 F 5 +2s 

3.583xl0' 9 

0.661 

11.06 

10.21 

From C2F4 

9 

s + C 2 F 4 —> 2s 

3.583xl0' 9 

0.661 

11.06 

10.21 

9,21 

10 

s + C 2 F 4 -> F + ^ 2 ^ + 2s 

3. 025x1 0‘ 9 

0.874 

16.41 

15.57 

9,21 

11 

s + C 2 F 4 . — > OF 2 + CF^ + 2s 

1.253xlO' 10 

1.514 

9.053 

14.33 

21 

12 

e + C 2 F 4 — > OF 3 + CF + + 2s 

5. 874x1 O' 9 

0.619 

19.29 

13.72 

9,21 

13 

s + O 2 F 3 — > C 2 F 3 " ^ 

3.583xl0' 9 

0.661 

11.06 

10.21 

From C2F4 

14 

e + CF 3 ->CF 3 +2e 

1.356xl0' 9 

0.796 

9.06 

10.0 

5 

15 

e + CF 3 — > F 4 -CF 2 + 2s 

7.020xl0 ' 9 

0.430 

16.28 

17.1 

5 

16 

S + CF 3 — > 2F + CF + + 2s 

4.148x10'® 

-0.341 

24.28 

21.4 

5 

17 

S + CF 2 — > CF 2 + 2 s 

1.103x10'® 

0.393 

11.37 

10.0 

5 

18 

e + CF 2 -> F + CF + + 2e 

5.434xl0 ' 9 

0.561 

14.29 

14.2 

5 

19 

s + CF->*CF + +2s 

5.480xl0' 9 

0.556 

9.723 

9.1 

5 

20 

e + F ->• F* + 2e 

2.397xl0 ' 9 

0.848 

16.78 

17.4 

9 

21 

s + F 2 - 4 - FJ + 2 s 

2 . 886 xl 0' 9 

0.881 

15.91 

5.7 

9 

22 

e+C ->C + +2s 

1.272xl0' 9 

1.648 

7.853 

11.3 

30 

Electron Impact Dissociation: 
23 e+C 4 F 8 ->2C 2 F 4 +e 

8.709x10-® 

0.042 

8.572 

2.16 

See text 

24 

S + C 4 F 3 — > CgFg + CF 2 + 6 

8.709xl0' 9 

0.042 

8.572 

3.25 

See text 

25 

s + CgFg — > C 2 F 4 + CF 2 + s 

1.070x10® 

0.230 

7.451 

4.53 

Avg. c-C 4 F 8 

26 

6 + C 3 F 5 — > C 2 F 3 + CF 2 + 6 

1.070x10'® 

0.230 

7.451 

3.16 

c 2 f 4 

Avg. c-C 4 F 8 

27 

e + C 2 F 4 — > 2 CF 2 + s 

1.315x10'® 

0.412 

6.329 

2.88 

c 2 f 4 

9 

28 

s + C 2 F 3 -> CF 2 + CF + s 

3.297x10'® 

0.412 

6.329 

2.73 

From C 2 F 4 

29 

s + OF 3 — > OF 2 + F + s 

6.484x10'® 

-0.959 

11.25 

3.8 

9 

30 

e + CF 3 — ► CF 2 + F + s 

7.941x10'® 

-0.452 

12.10 

3.8 

9 

31 

s + CF 2 — > CF + F 4 6 

8.113xl0' 9 

0.386 

8.739 

5.4 

See Sec. V 


32 

e + CF 2 -+C+2F + e 

1.390x1 0' 8 

-1.164 

49.87 

11.0 

6 

33 

s + OF — >C+F+e 

1.634x1 O' 8 

-0.002 

13.05 

5.6 

6, From CO 

34 

e + F 2 -> 2F+e 

1.080x1 0‘ 8 

-0.296 

4.464 

1.71 

9 

Dissociative Attachment: 






35 

© 4 - C 4 F 8 — > (^Fg CF + 

2.789x1 O’ 9 

-1.277 

5.392 

0.0 

9 

36 

e + CF 3 CF 2 + F" 

3.000xl0' n 

0.0 

0.0 

0.0 

See text 

37 

6 4- CF 2 — ^ OF + F 

3.000xl0" u 

0.0 

0.0 

0.0 

See text 

38 

© 4 - OF — >C + F 

3.000xl0' u 

0.0 

0.0 

0.0 

See text 

39 

6 4 - F 2 — > F + F 

5.768xl0' 9 

-1.465 

0.539 

0.0 

9 

Ion-Molecule Reaction: 






40 

c 4 f 8 +cf^c 2 f 4 +f+c 3 f 5 + 

2.100xl0' n 

0.0 

0.0 

0.0 

9 

41 

cf 3 +cf 2 + ->cf 2 +cf 3 + 

1.480x1 O' 9 

0.0 

0.0 

0.0 

9 

42 

cf+cf 2 -»c+cf 3 

2.060x1 O' 9 

0.0 

0.0 

0.0 

9 

43 

c+cf 2 ->cf+cf + 

1.040x1 O' 9 

0.0 

0.0 

0.0 

9 

44 

cf 3 +cf + ->cf+cfj 

1.710xl0' 9 

0.0 

0.0 

0.0 

9 

45 

cf 3 +c + -»cf+cf 2 + 

2.480x1 O' 9 

0.0 

0.0 

0.0 

9 

46 

CF + C + ->C + CF + 

3. 180x1 O' 9 

0.0 

0.0 

0.0 

9 

47 

cf 3 +f + ->F 2 +CF 2 + 

2.090xl0‘ 9 

0.0 

0.0 

0.0 

9 

48 

cf 2 +f + -»-f 2 +cf + 

2.280xl0‘ 9 

0.0 

0.0 

0.0 

9 

49 

cf+f + ->f 2 +c + 

2.7 10x1 0‘ 9 

0.0 

0.0 

0.0 

9 

50 

C + F + -> F + C + 

1.170xl0‘ 9 

0.0 

0.0 

0.0 

9 

Ion-Ion Recombination: 






51 

C 2 F 4 + 4 -F"->CF 2 4 -CF + F 2 

8.200xl0' 8 

0.0 

0.0 

0.0 

9 

52 

C3F5 +F — * C 2 F 4 +CF 2 

8. 000x1 O’ 8 

0.0 

0.0 

0.0 

9 

53 

CF 3 + + F'->CF 2 + F 2 

8.700xl0' 8 

0.0 

0.0 

0.0 

9 

54 

cf 2 + +f~ ->cf+f 2 

9.100xl0' 8 

0.0 

0.0 

0.0 

9 

55 

CF + +F~ -»CF+F 

9.800xl0‘ 8 

0.0 

0.0 

0.0 

9 

56 

C + +F'->-C+F 

1.200x1 O’ 7 

0.0 

0.0 

0.0 

9 

57 

F + +F~ ->F+F 

l.lOOxlO' 7 

0.0 

0.0 

0.0 

9 

58 

F 2 + +F‘-*F 2 +F 

9.400xl0‘ 8 

0.0 

0.0 

0.0 

9 

Detachment: 






59 

e 4- F -» F 4- 2e 

1.6153xl0' 7 

-0.436 

12.87 

3.0 

6 

Electronic Excitation: 






60 

6 4- C — ► C^D) 4- © 

4.357x1 O' 8 

-0.504 

2.12 

1.264 

31 

61 

e + C-» C( 1 S)4-e 

6.199xl0' 9 

-0.504 

3.59 

2.68 

31 

62 

e 4- F F( 4 P) + e 

3.544x1 O’ 9 

-0.650 

14.33 

12.70 

31 

63 

e + F-> F( 2 P)4-e 

3.618xl0' 9 

-0.053 

13.83 

12.99 

31 

Surface Recombination: 






64 

C+F->CF 

Recomb. Coeff 2y/(2-y)=0.2 



See Sec. V 

65 

cf+f->cf 2 

Recomb. Coeff. 2y/(2-y)=0.2 



See Sec. V 


